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The similarity renormalization group (SRG) is based on unitary transformations that suppress 
oft-diagonal matrix elements, forcing the hamiltonian towards a band-diagonal form. A simple 
SRG transformation applied to nucleon-nucleon interactions leads to greatly improved convergence 
properties while preserving observables, and provides a method to consistently evolve many-body 
potentials and other operators. 
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Progress on the nuclear many-body problem has been 
hindered for decades because nucleon-nucleon (NN) po- 
tentials that reproduce elastic scattering phase shifts typ- 
ically exhibit strong short-range repulsion as well as a 
strong tensor force. This leads to strongly correlated 
many-body wave functions and highly nonperturbative 
few- and many-body systems. But recent work shows 
how a cutoff on relative momentum can be imposed 
and evolved to lower values using renormalization group 
(RG) methods, thus eliminating the troublesome high- 
momentum modes d The evolved NN potentials 
are energy-independent and preserve two-nucleon observ- 
ables for relative momenta up to the cutoff. Such po- 
tentials, known generically as Viow/cj are more pertur- 
bative and generate much less correlated wave functions 
0, i, i, i, IE : vastly simplifying the many-body prob- 
lem. However, a full RG evolution of essential few-body 
potentials has not yet been achieved. 

An alternative path to decoupling high-momentum 
from low-momentum physics is the similarity renormal- 
ization group (SRG), which is based on unitary transfor- 
mations that suppress off-diagonal matrix elements, driv- 
ing the hamiltonian towards a band-diagonal form [^, 
[Ulill- The SRG potentials are automatically energy in- 
dependent and have the feature that high-energy phase 
shifts (and other high-energy NN observables), while typ- 
ically highly model dependent, are preserved, unlike the 
case with Viowk as usually implemented. Most impor- 
tant, the same transformations renormalize all operators, 
including many-body operators, and the class of trans- 
formations can be tailored for effectiveness in particular 
problems. 

Here we make the first exploration of SRG for nucleon- 
nucleon interactions, using a particularly simple choice of 
SRG transformation, which nevertheless works exceed- 
ingly well. We find the same benefits of Mowfc: more 
perturbative interactions and lessened correlations, with 
improved convergence in few- and many-body calcula- 
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tions. The success of the SRG combined with advances 
in chiral effective field theory (EFT) [13, opens the 
door to the consistent construction and RG evolution of 
many-body potentials and other operators. 

The similarity RG approach was developed indepen- 
dently by Glazek and Wilson 8] and by Wegner Q. We 
follow Wegner's formulation in terms of a flow equation 
for the hamiltonian. The initial hamiltonian in the center 
of mass H = T^d + V , where Trci is the relative kinetic 
energy, is transformed by the unitary operator U{s) ac- 
cording to 



H, = U{s)HU\s) = T, 



rcl 



(1) 



where s is the flow parameter. This also defines the 
evolved potential Vg, with Tj-ci taken to be independent 
of s. Then Hs evolves according to 



ds 



with 



^Gs) = ^(7t(s) = -,;t(,) 
ds 



(2) 



(3) 



Choosing r/(s) specifies the transformation. Here we 
make perhaps the simplest choice |To| . 



which gives the flow equation, 
dH, 



ds 



(4) 



(5) 



Other choices will be studied elsewhere p^ . 

For any given partial wave in the space of relative mo- 
mentum NN states, Eq. ([5]) means that the potential in 
momentum space evolves as (with normalization so that 
1 = -I \q)q'^ dq{q\ and in units where hj^/M = 1), 



dVs{k,k') 
ds 



^(fc2 _ k'^ fVs{k,k') 

f- Hq^dq (k^ + k'^-2q^) 
I" Jo 

xVs{k,q)Vs{q,k') . (6) 
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FIG. 1: (Color online) Contour plots of Vs{k,k') illustrating the evolution with A = s"^/* for ^So (left) and ^Si (right). The 
initial potential on the left is a chiral N^LO potential with a 600 Me'V cutoff and on the right is an N^LO potential with a 
550 MeV cutoff on the Lippmann-Schwinger equation and a 600 Me'V cutoff on a regularized spectral representation of two-pion 
exchange [l3( |. 
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FIG. 2: (Color online) S-wave phase shifts from the two chiral EFT N'^LO potentials from Fig. [T] For each initial potential, 
the phase shifts for different A agree to within the widths of the lines at all energies shown. 



(The additional matrix structure of Vg in coupled chan- 
nels such as '^Si~^Di is implicit.) For matrix elements 
far from the diagonal, the first term on the right side of 
Eq. ([6]) evidently dominates and exponentially suppresses 
these elements as s increases. The parameter A = s^^/"* 
provides a measure of the spread of off-diagonal strength. 
While further analytic analysis is possible, we instead 
turn to a numerical demonstration that the flow toward 
the diagonal is a general result. By discretizing the rel- 
ative momentum space on a grid of gaussian integration 
points, we obtain a simple (but nonlinear) system of first- 
order coupled differential equations, with the boundary 
condition that Vs{k, k') at the initial s (or A) is equal to 
the initial potential. 

The evolution of the hamiltonian according to Eq. ^ 
as s increases (or A decreases) is illustrated in Fig. (TJ 
using two initial chiral EFT potentials [l^, • On the 
left in Fig. [1] is "'^Sq starting from the harder (600 MeV 
cutoff) potential from Ref. [l^l, which has significant 
strength near the high-momentum diagonal, and on the 



right is the S-wave part of the '^Si-^Di coupled channel 
starting from one of the potentials from Ref. [3], which 
has more far off-diagonal strength initially and compar- 
atively weaker higher-momentum strength on the diago- 
nal. The initial momentum-space potential differs signif- 
icantly among interactions that are phase-equivalent up 
to the NN inelastic threshold, but these examples show 
characteristic features of the evolution in A. In particular, 
we see a systematic suppression of off-diagonal strength, 
as anticipated, with the width of the diagonal scaling as 
A^. The same behavior is observed when evolving from 
conventional high-precision NN potentials, such as Ar- 
gonne uis , or other (softer) chiral potentials [3, EBl • 

Since the SRG transformation is unitary, observables 
are unchanged at all energies, up to numerical errors. 
This is shown by Fig. [21 in which phase shifts for the two 
chiral EFT potentials are plotted, including the values 
at high energies where they are not constrained by data 
(above E'lab = 300 MeV). For a given potential, there is 
no visible variation with A. Similarly, the binding en- 



3 



ergy and asymptotic normalizations for the deuteron are 
independent of A [15]. 

As A is lowered, different initial potentials flow to sim- 
ilar forms at low momentum while remaining distinct at 
higher momentum. The low-momentum parts also be- 
come similar to Viowfc potentials. These observations are 
illustrated in Fig. [3] for two particular slices of the poten- 
tials from Fig. ^ They will be explored in much greater 
detail in Ref. 11411 . 
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FIG. 3: (Color online) Matrix elements of the evolved SRG 
potentials at A = 2 fm~^ for ^So (top, diagonal elements) 
and ^Si (bottom, off-diagonal elements) for the same initial 
potentials as in Fig. [T] Also shown is the Viowfc potential 
with a smooth (exponential) regulator for momentum cutoff 
A = 2 fm~^, evolved from the two potentials (600 MeV above 
and 550/600 MeV below). 

We can quantify the perturbativeness of the potential 
as we evolve to lower A by using the eigenvalue analysis 
introduced long ago by Weinberg [l^ and recently ap- 
plied in an analysis of Mow k potentials [B] . Consider the 
operator Born series for the T-matrix at energy E (for 
simplicity we assume -E < 0): 



TiE)^Vs + V,- 



1 



-Vs + 



E — Trcl 

By finding the eigenvalues and eigenvectors of 
1 



E — Trcl 

and then acting with T(E) on the eigenvectors, 

r(i?)|F,) = IF.) (1 + 77^ + 772 + ...) , 



(7) 



(8) 



(9) 



it follows that nonperturbative behavior at energy E is 
signaled by one or more eigenvalues with jry. (_E)| > 1 
(See Ref. Q for a more detailed discussion in the 
context of evolving Viowfc potentials.) 

It suffices for our purposes to consider a single energy 
(e.g., E = 0), and to consider only the negative eigenval- 
ues, which are associated with the short-range repulsion. 




FIG. 4: (Color online) The largest repulsive Weinberg eigen- 
values as a function of A in the ^So channel and the ^Si-^Di 
coupled channel for the same initial potentials as in Fig. [T] 



Weinberg eigenvalues at zero energy are shown as a func- 
tion of A in Fig. Hfor the ^So channel and the '^Si-'^Di 
coupled channel. In both channels, the large negative 
eigenvalues at large A refiect the repulsive core of the ini- 
tial potentials. They rapidly evolve to small values as A 
decreases to 2fm~^ and below, as also observed with the 
corresponding Viowk evolution 0]. However, the interme- 
diate increase for the sub-leading eigenvalues in -'^Sp is a 
new feature of the SRG that merits further study. 
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FIG. 5: (Color online) The absolute error vs. A of the pre- 
dicted deuteron binding energy from a variational calculation 
in a fixed-size basis of harmonic oscillators (A'^max^t^ excita- 
tions). The initial potential is from Ref. [131] . 

The more perturbative potentials at lower A induce 
weaker short-range correlations in few- and many-body 
wave functions, which in turn leads to greatly improved 
convergence in variational calculations. This is illus- 
trated via calculations of the binding energy of the 
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deuteron and triton by diagonalization in a harmonic os- 
cillator basis, as shown in Figs. [S] and [SI For a fixed basis 
size, a more accurate estimate is obtained with smaller A 
or, conversely, at fixed A the convergence with basis size 
becomes more rapid. The improvement in convergence 
is similar to that found with smoothly regulated Mow fc 
potentials evolved fi'om chiral N'^LO potentials frj. At 
finite density, analogous effects led to perturbative be- 
havior in nuclear matter for Viow k potentials [1| . Results 
for the G matrix support a similar conclusion for the SRG 
potentials (see the website in Ref. Q for pictures). 




FIG. 6: (Color online) The variational binding energy for 
selected A of the triton with two-nucleon interactions only, as 
a function of the size of the harmonic oscillator space {Nmaxhuj 
excitations), for the same initial potentials as in Fig. [1] 

In Fig. [6l the calculations for different A converge to 
different values for the binding energy of the triton. This 
reflects the contributions of the omitted (and evolving) 
three-body interaction. The evolution with A of the bind- 
ing energy with NN interactions only, which is also the 
evolution of the net three-body contribution, follows a 
similar pattern to that seen with Mow k 0, 0] : a slow de- 
crease as A decreases, reaching a minimum for A between 
1.5 fm~^ and 2 fm~^, and then a rapid increase. 

The consistent RG evolution of few-body interactions 
is an important unsolved problem for low-momentum po- 
tentials. In Mow k calculations to date, an approximate 
evolution is made by fitting the leading chiral EFT three- 
body force at each cutoff while evolving the two-body 
interaction exactly Q. Generalizing the RG evolution 
to include the three-body interaction is, at least, tech- 
nically challenging. This is because the machinery used 
to construct Mow k requires the solution of the full three- 
nucleon problem (i.e., bound state wave functions plus 
all scattering wave functions in all breakup channels) to 
consistently evolve the three-nucleon interactions. 

In contrast, the SRG evolution follows by directly ap- 
plying Eq. (O in the three-particle space. To show the 
basic idea, we adopt a notation in which V12 means the 



two-body interaction between particles 1 and 2 while M123 
is the irreducible three-body potential. In standard treat- 
ments of the nuclear three-body problem (for example, 
Ref. [l3l)j the relative kinetic energy for equal mass par- 
ticles with mass m is decomposed as (with the total mo- 
mentum K = ^ . ki = 0): 



rcl 



f ^ 9ttj IT) Am 



Z = 1,2,3, (10) 



where I denotes which set of Jacobi coordinates are being 
used to describe the relative motion (i.e., I = 1 means 
Pi \{^2 - ka) and qi = |[ki - i(k2 + ka)], and so 
on). In the notation here. 



Hi 

m 



t23 



and 



3q? 

4m 



Ti , 



(11) 



and similarly for I — 2 and I — Z. 

Now we start with the hamiltonian including up to 
three-body interactions (keeping in mind that higher- 
body interactions will be induced as we evolve in s but 
will not contribute to three-body systems): 



Hs = T,el + M12 + Mi3 + M23 + M 



23 



T 



rcl 



Ms 



(12) 



(Note: all of the potentials depend implicitly on s.) The 
relative kinetic energy operator Trci can be decomposed 
in three ways: 

Tlcl = T12 + T3 = Ti3 + T2 = r23 + Ti , (13) 

and Ti commutes with Vjk, 

[T3, M12] = [r2, M13] = [Ti, M23] = , (14) 

so the commutators of Trci with the two-body potentials 
become 



[T;el,Mi2] = [Ti2,Mi2] 



(15) 



and similarly for M13 and M23. 

Since we define Trei to be independent of s, the SRG 
equation ([5]) for the three-body Hamiltonian Hs simpli- 
fies to 



dM 
ds 



dVi2 , dMi3 , dM23 , dVi 



23 



ds 



ds 



ds 



ds 



= [[T,.el,M],ff, 



(16) 

with M, from Eq. p2p . The corresponding equations for 
each of the two-body potentials (which are completely 
determined by their evolved matrix elements in the two- 
body systems, e.g., on a gaussian momentum grid) are 



dV^ 



12 



ds 



= [[ri2,Mi2],(Ti2+Mi2)] 



(17) 



and similarly for M13 and M23. After expanding Eq. ([Tf 
using Eq. (fT^ and the different decompositions of T^ci: 
it is straightforward to show using the equations for the 



5 



two-body potentials that the derivatives of two-body po- 
tentials on the left side cancel precisely with terms on 
the right side, leaving 



dV-i 



123 



ds 



[[^12, 1^12], (T3 + Vi3 + V23 + V123)] 

+ {[Ti3, V13], (T2 + Vu + V23 + 1/123)] 
+ [[T23, V23], (Ti + V12 + Vu + V123)] 
+ m,uV,23],H,] . (18) 



The importance of these cancellations is that they elim- 
inate the "dangerous" delta functions, which make set- 
ting up the i nteg ral equations for the three-body system 



problematic [17[. We emphasize that the s-dependence 



of the two-body interactions on the right side of Eg.fTSlis 
completely specified by solving the two-body problem in 
Eq. ^ . This is in contrast to RG methods that run a cut- 
off on the total energy of the basis states (e.g., Lee-Suzuki 
and Bloch-Horowitz techniques). Such methods generate 
"multi-valued" two-body interactions, in the sense that 
the RG evolution of two-body operators in ^ > 2 sys- 
tems depends non-locally on the excitation energies of 
the unlinked spectator particles [H, [l^ . 

Further simplications of Eq. (fT5|) follow from sym- 
mctrization and applying the Jacobi identity, but this 
form is sufficient to make clear that there are no discon- 
nected pieces. The problem is thus reduced to the techni- 
cal implementation of a momentum-space decomposition 
analogous to Eq. Note in particular that the evolu- 
tion of V123 is carried out without ever having to solve 



a bound-state or scattering problem. We have verified 
that this formalism does generate three-body interactions 
that leave eigenvalues invariant for simple model hamil- 
tonians, such as a two-level system of bosons. Work is in 
progress on proof-of-principle tests using one-dimensional 
many-body systems (to avoid angular momentum com- 
plications) and the three-dimensional nuclear problem 
will be tackled soon. 

In summary, the SRG applied to nuclcon-nucleon po- 
tentials works as advertised even for a simple choice of 
transformation, driving the hamiltonian (in momentum 
space) towards the diagonal, making it more perturbative 
and more convergent in few-body calculations. There is 
much to explore, such as the nature of the decoupling 
of high- and low-energy physics implied by Fig. [T] (see 
Ref. ^5|] for demonstrations of decoupling and the evolu- 
tion of an operator in the two-body space) and whether 
other choices of 77(3) instead of Eq. ([4]) could be more 
effective in making the hamiltonian diagonal. For exam- 
ple, the replacement Tici Hd, where Hd is the diagonal 
part of the hamiltonian, or some function of Tj-ci are easily 
implemented. Most important is the consistent evolution 
of nuclear three-body operators [T3 |. 
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